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We investigate dynamical many-body localization and delocalization in an integrable system of 
periodically-kicked, interacting linear rotors. The linear-in-momentum Hamiltonian makes the Flo- 
quet evolution operator analytically tractable for arbitrary interactions. One of the hallmarks of 
this model is that depending on certain parameters, it manifests both localization and delocaliza¬ 
tion in momentum space. We present a set of “emergent” integrals of motion, which can serve as 
a fundamental diagnostic of dynamical localization in the interacting case. We also propose an ex¬ 
perimental scheme, involving voltage-biased Josephson junctions, to realize such many-body kicked 
models. 


I. INTRODUCTION 

Recently, there has been a lot of interest and progress 
in understanding Anderson-type localization properties 
of disordered, interacting many-body systems. Notably, 
the remarkable phenomenon of many-body localization 
(MBL) was discovered 1-8 . In an isolated system, MBL 
manifests itself in the localization of all eigenstates and 
leads to the breakdown of ergodicity and violation of 
the eigenvalue thermalization hypothesis 9,10 , forcing to 
revisit the very foundations of quantum statistical me¬ 
chanics 11,12 . 

In this work, we ask whether a driven interacting sys¬ 
tem can be dynamically many-body localized. We answer 
this question in the affirmative and present an exactly 
solvable model of a kicked chain of interacting linear ro¬ 
tors, which shows both dynamical MBL and delocalized 
regimes. 

A quantum kicked-rotor is a canonical model of quan¬ 
tum chaos 13,14 which exhibits dynamical localization in 
momentum space. The time dependent Schrodinger 
equation for a general kicked rotor is given by (here and 
below, we set h = 1 and the driving period T = 1), 

i d t ip(6,t) = \2Tta{-\d s ) 1 + K(6)6(t - n)]ip(0,t). (1) 

In the ground breaking paper 13 , Fishman, Prange, and 
Grempel proved that the eigenvalue problem for the 
Floquet operator of a kicked rotor is equivalent to 
that of a particle hopping in a (quasi)periodic potential 
((ir)rational a) in one-dimension 1 ' 3 given by, 

W m + r u r + tan [oj — 2iiam l ]u m = Eu m . (2) 

r 

This mapping traced the origin of the dynamical localiza¬ 
tion in driven systems to Anderson localization in time- 
independent settings. While the quadratic rotor (l = 2) 
is nonintegrable, the linear rotor model (l = 1) in Eq. (1) 
was exactly solved in Refs. 13,15 “ 18 . The corresponding 


integrable lattice version in Eq. (2) with l = 1 is dubbed 
Maryland model (MM) 19,20 (See section Al for techni¬ 
cal details of this mapping). For the linear rotor, both 
classical and quantum dynamics is integrable, and the 
dynamical localization (absence of chaos) is due to the 
existence of a complete set of integrals of motion 18 . How¬ 
ever, the incommensurate MM is an Anderson insulator 
with no classical interpretation 1 ' 1 even though the linear 
rotor manifests classical integrability. Thus, the dynam¬ 
ical localization for the quantum version of both linear 
(l = 1) and quadratic rotor (l = 2) seems to stem from 
the Anderson mechanism 15 . 

In this work, we generalize the linear kicked rotor 
model by considering an interacting chain of driven rotors 
in order to understand the anatomy of many-body local¬ 
ization in the dynamical space. One of the remarkable 
features of this many-particle model is that it manifests 
both localization and delocalization in dynamical space 
depending on whether the components of the a are irra¬ 
tional or rational. 

Let us emphasize from the outset that the model we 
consider is integrable for all parameters owing to the 
first order differential operator. This leads to integrals 
of motion (IOMs), which are local in the spatial (angu¬ 
lar) variables. The underlying integrability is a special, 
non-universal feature of our model (3), which allows us to 
solve it exactly. However, the existence of the local-in-0 
IOMs has no direct relation to dynamical MBL, which 
occurs in momentum space. 

As shown below, dynamical MBL is accompanied by 
the appearance of additional integrals of motion bounded 
in momentum space - a central result of this work. The 
non-interacting version of these additional IOM’s for the 
linear kicked rotor was pointed out by Berry in Ref. [ 18]. 
As argued below, these “emergent” IOMs and the dy¬ 
namical MBL are universal phenomena, which would sur¬ 
vive in non-integrable generalizations of the model (which 
however is not analytically solvable). 

To capture the dynamical localization for this interact- 
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ing model we monitor three indicators: energy growth at 
long times, (momentum degrees of freedom at long times, 
and the existence of integrals of motion in the momentum 
space. Below, we first present the detail of our model in 
Section II. The main analytical results are outlined in 
Section III, and their numerical analysis and key con¬ 
clusions are given in Section IV. An outline of technical 
derivation of the results is given in Section V.The exper¬ 
imental proposal to realize our model (3) is explained in 
Section VI. In section VII we list our conclusions and 
Section VIII contains acknowledgements. Finally, Ap¬ 
pendix A1 and Appendix A 2 are devoted to the sum¬ 
mary of the connection of the rotor problem to a lattice 
model and the correspondence between our model and a 
d-dimensional lattice, respectively. 

Throughout the text, in all the summations JV ... or 
"Yhij ■■■ we only consider i ^ j. So that, expressions like 
JV 1 /(ai — ctj ) are not divergent. The denominator van¬ 
ishes only when there is a resonance, such as a* —> aj. 


II. THE MODEL 

We consider a many body interacting generalization, 

oo d 

H(t) =H 0 + V E d(t — n), with V = AT(d, ; ), 

n ——oo i =1 

d ^ d— 1 

H 0 = 27t ^2 + o E Jjj{Qj ~ Oj ) ( 3 ) 

i=l 

of the linear rotor model. 

Hq is the static Hamiltonian describing d particles on 
a ring, each rotating 27rai radians per one period of the 
kick. Oi is the position operator for the Ah particle on 
the ring and pi is its angular momentum operator, which 
has integer eigenvalues in the H = 1 units. These d parti¬ 
cles interact through a translationally invariant two-body 
potential Jij{0i — Oj). This form of the interaction en¬ 
sures conservation of momentum . The rotors are driven 
by V(t) which contains periodic delta function impulses, 
where the strength of the impulse is given by a general 
periodic one-body potential K(0i). The local potentials 
are periodic with 2n and therefore the most general form 
can be written as K{0j) = ■ Here, k m is the 

mth Fourier component of the potential that acts on the 
jth particle. The periodic form of the interaction poten¬ 
tial J can be written as J\j = Yh m b l J l e lm ( ei ~ e U. Here, 
b% is the mth Fourier component of the interaction po¬ 
tential between the Ah and jth particle. Reversing i and 
j and replacing m with —m in this Fourier expansion has 
no effect on the components, therefore 6^ = bU Tn . This 
property will be handy while deriving formulas through¬ 
out the text. 

Note that in our model the localization is a conse¬ 
quence of incommensurate driving period and angular 
velocities a. This situation is different from recent works 


using Floquet analysis to probe dynamical properties of 
MBL states of disordered Hamiltonians 2 1-24 (in these pa¬ 
pers, a Floquet perturbation is imposed on a state, many- 
body localized in coordinate space, while our goal is to 
induce dynamical many-body localization in momentum 
space by the Floquet perturbation). 


III. MAIN RESULTS 


A. Energy dynamics 


Following the conjecture of D’Alessio and 
Polkovnikov 25 , we test the dynamical localization 
by computing the energy growth as a function of time at 
long times. The average energy after N kicks (equivalent 
to time) can be written as, E(N) = (;i^n\Ho\Un) ■ To 
compute it, we write 

I Un) = up IV’o), Uf = e-'V^ 0 , (4) 


where Uf is the Floquet evolution operator, which cap¬ 
tures the state of the system after each kick. Owing to 
the linear dependence of the Hamiltonian on the momen¬ 
tum, we can explicitly compute this expectation. 


d 

E(N) = E( 0) + EE 27ra,;(f m j)o 

2 = 1 m 


sin^mNirai) 

sin(mnai) 


( 5 ) 


In the above expression, E( 0) = (^ol^ol^o) corresponds 
to the average many body energy over the initial state. 

r mi = —i mk rn e lm ^ i+7TO ‘ i ^ N+1 ^ depends on the chosen 
form of K{6) averaged over the initial state and due to 
its periodic nature is a bounded function of the number 
of kicks N. The growth of energy for large N is then 
completely dependent on the nature of ai appearing in 
the ratio sln (( riJV7rQ: ») Note that this expression is com- 
pletely independent of interactions, which we prove in 
Section V A. This is a consequence of momentum con¬ 
serving interactions and this property is no longer valid 
when the translational invariance of interactions is bro¬ 
ken. 


B. Momentum dynamics 

Another indicator for dynamical localization is the 
spread in the momentum degrees of freedom. The Ah 
momentum after N kicks, Pi(N) = (UnIPiIUn) is given 
by the following expression, 


Pi{N) 


<Pi)o+E< f mi) 0 

m 


sm(m.Nnai) 

sin(m7rcq) 


+E< f ™> 


sm(mNnAaij) 

mnAaij 


( 6 ) 
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We have defined Aaij = on — ay. In the above expression, note that there are cases where moti can get arbitrar- 

the first term is the ith momentum in the initial eigen- ily close to integers (Liouville numbers) and total energy 

state (pi )o = (4’o\Pi\' t Po)■ The second term corresponds and momentum are no longer bounded. In the single ro- 

to the kicking potential as defined in Eq. (5). The last tor case, such a situation yields interesting consequences 

term depends explicitly on the form of interaction via like marginal resonance and a mobility edge in the mo- 

f“b = -i rn \ J i3 i ^m(e i -e :j +-KNAa ij ) and is a bounded func- mentum lattice, 18 ’ 31 . In this study we exclude this pos¬ 
tion of N. The growth of momenta at long times cor- sibility and refer to generic irrationals only, 

responding to the last term is completely determined by 
the ratio sin( ™j^ a < 

A. Case I: qi ... otd are distinct generic irrationals 


C. Integrals of motion 

Recent works have shown that the existence of inte¬ 
grals of motion (IOM) can be used as a diagnostic to 
quantify both non-interacting Anderson localization 26,27 
and many body localization 28-30 . In the context of dy¬ 
namical localization for the model at hand, we work in 
the momentum basis and search for existence of IOMs in 
this basis. We begin by constructing IOMs by identifying 
operators C* that satisfies [C*, Uf\ = 0 and [C,;, Cj] = 0. 
The existence of these IOMs encode information about 
dynamical localization. Operators that commute with 
Uf satisfy the property (C'»)jv+i = ( Ci)w and are given 

by 

1 rnfc m e im ^ +7ra b 1 

1 2 sin(m7rai) 2 n(ai — a,-) 

m v 7 mj J 

( 7 ) 

This expression is a generalization of the constant of mo¬ 
tion given by Berry 18 . The derivation of this expression 
is given in Sec. VB. Since C* is an IOM, pi is bounded in 
time as long as the series in the last two terms converges. 
The delocalization of pi with time will occur due to the 
diverging denominators in the 9 dependent terms. 

We reiterate that the model always contains d inte¬ 
grals of motion, Bi = 0i/cti (mod 27r) {[Bi,Up] = 0) 
which results in integrability for both localized and de¬ 
localized cases. For example, if a, = r.i/si is rational, 
the existence of Bi results in integrability even though 
C\...d do not exist. Since Bi is momentum independent, 
its existence cannot bind pi(N). For the fully localized 
case, we have additional d IOM’s C\...d given in Eq. (7) 
that explicitly depend on p t and thus constrain it. Thus, 
Ci s can be thought of as “emergent” IOM’s constraining 
the momentum growth, resulting in dynamical MBL and 
“emergent” superintegrability. 

IV. DYNAMICAL MANY-BODY 
LOCALIZATION AND THE STRUCTURE OF a 

Equipped with three analytical expression for the en¬ 
ergy growth, momentum growth and IOMs [ eqs. (5) 
to (7)], we now diagnose the dynamical localization for 
three distinct cases with varying structure of a. We 


For irrational values of a*, the total energy in Eq. (5) 
is always a bounded function of N since moq ^ integer. 
In Fig. 1(AI), we fix initial states to be momentum eigen¬ 
states ipQ (6) which is e lp °' 6 up to a normalization factor. 
For momentum eigenstates, (T m i)o = 0 as the expec¬ 
tation value involves the integral f dde lm0 over a circle. 
Thus, we plot the root mean square (RMS) of system’s 

energy a(E) = y/(H$) - ( H 0 ) 2 as a function of N. Fig. 1 
(AI) shows boundedness in the spread of the energy as 
a function N. The momentum growth shown in Eq. (6) 
has contributions from both interactions and the kick¬ 
ing potential. In Fig. 1 (All) we plot RMS deviation of 
momenta cr(pi)... a(pd) (we used d = 10 for the specific 
simulation) and show that the spread in the momenta is 
bounded as a function of N. For this case the integrals 
of motion in Eq. (7) exist and convergent. Thus, all the 
diagnostics for this case point towards a true many body 
dynamical localization. 


B. Case II: a i = «2 and a 2 ■ • ■ Ud are distinct generic 
irrationals 

For this case, the energy remains a bounded func¬ 
tion of N since moq are not integers. Thus, the RMS 
deviation cr(E) is bounded as shown in Fig. (1 BI). 
However, the second term in the momentum expres¬ 
sion (Eq. (6)) develops a resonance (for the momenta 
Pi and P2) since oq —► a 2 . Due to this resonance 
term, sin(mV7rAai2)/(m7rAai2) ~ N as oq —► oq. 
This resonant growth is reflected in the RMS deviation 
of cr(pi) and er(p 2 ) growing linearly with N , while the 
■ ■ ■ <?{pio) remain bounded as shown in Fig. 1 (BII). 
This is a striking scenario where the resonant momenta 
are not localized even if the total energy is bounded for 
large N. However, the delocalization of p\ and P2 in time 
is reflected in the break down of IOMs C-\ and C 2 due to 
diverging denominators in the resonant limit oq —> a2- 
For this case, the bounded total energy fails to diagnose 
delocalization as shown in Fig. 1 (BI, BII). We note that 
this scenario has no analogue in the non-interacting limit. 
Notice that interactions we considered possess transla¬ 
tional invariance, i.e. in the form J{9i — 9j) given in 
Eq. (3) and therefore interactions conserve momentum. 
The dichotomy between the energy growth and momen¬ 
tum growth is a result of conservation of momentum and 
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(AI) a £ irrational, cti ctj (BI) a £ irrational, ct\ = ai (Cl) «i = 1/2, 02,..,10 £ irrational 
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FIG. 1. (Top panel) Plots showing evolution of the root mean square deviation of energy, a(E) = yj(H q) — {Ho} 2 , as a function 
of time, labeled by the number of kicks, N. (Bottom Panel) Plots showing evolution of the root mean square deviation of 
individual momenta in an ensemble of 10 particles cr(pi) = yj (p?) — (Pi ) 2 as a function of time, labeled by the number of kicks, 
N. The initial states of the rotors are assumed to be definite momentum states. The total number of particles in this interacting 
ensemble is d = 10 and circular boundary conditions apply. We consider two body interaction term to be Jij = cos (8i — 9j). 
The periodic kicking potential is given by K(9) = J/m=i km cos(md) for all particles. The m-tli Fourier coefficient is fixed by 
km = 2 /m 2 , where z is a random complex number with real and imaginary parts uniformly distributed in the interval [0,1]. 
We make sure that K(9) is real. We choose a corresponding to three scenarios: Figs. (AI, All): If ip denotes the golden ratio, 
Qj = ( j/d)ip — (l/2)i p are all irrationals for j = 1... 10. Figs. (BI, BII): We consider Qj distinct irrationals as in Case A, but 
set cti = « 2 , which results in the resonant growth of cr(pi) and < 7 (^ 2 ), while cr(E) is bounded. Figs. (Cl, CII): We consider ctj 
as in Case A and set an = 1/2, which results in the resonant growth of &(E) and cr(pi). 


the linear dependence of energy on momentum. If we al¬ 
low interactions that break translational invariance, mo¬ 
mentum is no longer conserved and resonance due to in¬ 
teractions triggers unbounded growth (delocalization) of 
both momenta and the total energy. 


C. Case III: an = 1/2 and 012 ... a<j are distinct 
generic irrationals 

For rational ot\ = 1/2, the system develops a different 
kind of resonance compared to Case II. We consider a 
kicking potential with (fc a 7 ^ 0). The resonance condition 
inai £ integer can be satisfied for m = 2 and the ratio 
S sin(m 7 ra“)^ g rows as N. It results in energy growth and 
delocalization of momenta pi as shown in the RMS de¬ 
viations in Fig. 1 (Cl, CII). This is a converse situation 
to Case II where the energy delocalizes with the delocal¬ 
ization of pi even though the momenta P 2 ,P 3 , ■ ■ -Pio are 
bounded. This situation is again captured by the IOMs, 
where Ci does not exist, while C 2 ■ ■ ■ C 10 are well defined 
and convergent as seen in Eq. (7). 

We can generalize the above representative cases. For 


each rational a^, its integral of motion Ci breaks down 
and the corresponding momentum and energy diverge 
with time. For each pair oti = otj both C) and Cj break 
down and the corresponding momenta diverge in opposite 
directions, therefore Ci + Cj is still an IOM and the total 
momentum and energy of the pair are bounded. 

Other than the behaviour of energy in Case II in Sec¬ 
tion IV B, the rest of our analysis apply equally to the 
interactions that break translational symmetry. 


V. DERIVATION OF MAIN RESULTS 

Having established the physical understanding of local¬ 
ization for this model, we now sketch the brief derivation 
leading to the final results in eqs. (5) to (7). In the fol¬ 
lowing we consider the expectation of a generic operator 
as a function of time, X(N) = (i/jn\X\i/jn)■ We write 
the explicit evolution of the many body wave function 
between two successive kicks, |Vbv) = e~ lV e~ lH °\ipN-i)■ 
Notice that Hq = 27ra • p + \ Jij0i ~ 9j) con¬ 
tains the many-body interaction term which may seem 
daunting, however, the linear momentum term allows 
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a factorization in the Floquet operator. The Baker- 
Campbell-Hausdorff formula Z = ln(e A 'e Y ) becomes 
tractable when [X, Y] = sF. In this case, the result 
simply reads Z = X + sY/( 1 — e~ s ). Now fix m and 
let X = i27r(aipi + 02 ^ 2 ) and Y = iexp(im[$i — $ 2 ]) 
where Aa ±2 = a\ — ct 2 - The result of the commutator 
reads [X,Y] = i27rAai2mF, precisely the tractable case 
discussed above. Replacing s with i27rAai2m in the for¬ 
mula Z = X + sY/{ 1 — e~ s ) and inverting both sides of 
the equality, we obtain the following factorization 


i2nct-p— ib))) exp(im[#i —$ 2 ]) 


_ g-i&m exp(im[0i-S2]) e -i27ra-p (g( 


For this to hold, the Fourier coefficients must satisfy 


b ij = 


sin(7rmAa ii ) im7rAa ., 
irmAaij 171 


( 9 ) 


This argument can be generalized to more particles and 
Fourier coefficients. Due to linearity summations over 
particle indices i,j and Fourier indices m are introduced. 
All in all, we can factorize the evolution operator for our 
model in Eq. (3) as, 


\ip N ) = e 


— V-%J2a Jij(6i-Gj) p -i2ir3-p 


iVdv-r), (10) 


where we have defined the modified interaction term as, 


(11a) 

m 

= y hi, (nb) 

„ nrn(Aaij) 


The advantage of this factorization is that the operator 
e ^i27r a-p j s a translation operator in the position basis. 
We can rewrite Eq. (10) in the position basis by acting 
with (9 1 from the left. We define (0 |^>jv) = ^jv( 0) and 
express (9\e~ l27Ta 'P\ipN-i) = ipN-i{0 — 2nd). For a single 
kick we then have, 

ij> N {6) = e - iv{S) -'^ _ 27ra). (12) 

The above equation can be recursively iterated to yield, 


-i E [V(e-2-Kna)+^J2ij Jij(G-2irn3)} 

Wn(p) = e ™=° 


x ipo(9 ~ 2TrNa). (13) 


Here (9 — 27 rna) is a short hand notation for J l3 (0i — 

9j — 2irnAaij). 


A. Derivation of the evolution of energy and 
momentum averages 

Now consider a generic operator X = 
X(pi...pd',0i...9d)- The expectation value of this 


operator after N kicks is 

X(N) = J dO r N 0)X ( W d ’ e '- 9 *) ( 14 ) 

If we substitute X =pk, we get 
pi(N) = (pi) o- 

N / i _ _ 

'E ( diV(9 + 27 rna) + di~ J,j(0 + 2irna) 

n=\ \ 2 ij 

(15) 

Here, di = d/d9i and (...) 0 = / d0...|V^o| 2 - The contribu¬ 
tion from the kicking potential is: 



—di E V(9 + 27mcE) 



sin(?nA7raj) 

sin(?n7ra;) 

(16) 


The contribution from the interaction potential is 


~9i E b E 4 - 0j + 2™[ati - «j]) 


(17) 


= E^-)o 


sin(?nA7rAaji) 

nmAaii 


Putting together above expressions, we obtain the ex¬ 
pression for the momentum dynamics presented in 
Sec. (Ill B). 


Z -* G1 


E<^)« 


sin(m7raj) 

sin(miV7rAaij) 

mnAan 


(18) 


Now we derive the expression for the energy growth. 
Substituting A' = H 0 in Eq. (14), we have 

E(«)=£(o )+ ^2 m . ( f m .)o sill( ’; iWm ; ) 

sinirmraz) 

mi 

+ - E ](Jjj0i ~ + 2ttNA aij) — Jij(9i — 9j)) o 


E 2ra *< f ^) 0 


sin(mA7rAai :/ ) 

rmrAaa 


(19) 


In the above equation, a cancellation occurs between the 
interaction terms in the last two lines of Eq. (19) owing 
to the following relation, 

— E<4j$ — T 27rAAajj) — Jij(9i — 9j ))o 
ij 

/Wnt\ sin(mJV7rAa y ) 

--y j 2 TOi (r „ iJ ) 0 m . (20) 
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This completes the derivation of the energy growth 
shown in Sec. (19) 


d 

E(N) = E{ 0) + EE 27rai(f m i)o 

i —1 m 


sm(mNnai) 

sin(TO7raj) 


( 21 ) 


The dropping out of interaction terms from the total 
energy can be interpreted in the following way. As seen 
from Eq. (18), the contribution of interaction to the mo¬ 
mentum average picks up a negative sign when i and j are 
interchanged. This means interaction transfers momen¬ 
tum from one particle to the other at each kick, in other 
words momentum is conserved for each couple of rotors. 
Since the energy is linear in momenta, when the mo¬ 
menta of rotors are summed, the contribution of interac¬ 
tions vanishes. We emphasize that when the interactions 
break translational invariance, they no longer conserve 
momentum and in that case energy growth depends on 
interactions too. 


B. Construction of integrals of motion 


localization C\,,d are intact. In this case, of particular im¬ 
portance is the following combination of the integrals of 
motion 


H f = £>^<3,. ( 24) 

i 


Where H F is known as the Floquet Hamiltonian and is 
defined as, 


e ~ iIlF = = U F . (25) 

The quasienergy wavefunctions ‘ip u are simultaneous 
eigenstates of H F and CVs. The quasienergy-w state cen¬ 
tered around momenta (p) = M is 





fc m e im( ^ +7r “i ) 
2 sin(?n7ray) 


E 

ijm 


bw 




47 t m(a.i — ctj) 


(26) 


In this section, we outline the derivation involved in 
the construction of integrals of motion. By inspecting 
Eq. (15) and using the identity sin(m7ra) = [exp(im7ra) — 
exp(—im7ra)]/(2i), we can write 


pi(N + 1) — pi{N) = 

1 \ - mnh 


sin(m7TQf/) 


(J im(Gi+ivati)\ _ / e im(0i+nai) 
N 


N+l 


1 

2 i 7rA ai 


g i miOi-Oj)} _ ( e imiQi-Oj) 
N 


N+l 


( 22 ) 


noting that, the expression is valid whenever the denom¬ 
inators sin(m7raz) and A aij = ai — ay are non-zero, in 
other words, whenever the resonances are avoided. The 
above expression can be organized in a way that it man¬ 
ifests the IOMs, 

/ 1 y' rnk m e im ( §i+7rai) 1 \ 

\ 2 2—/ smimircti) 2 2—/ tHq- — ay) / 

\ rn mj J 

/ rnfc m e im ^ i+7rai ) 1 b^e irn ^-+ \ 

\ 1 2 2—/ sin(m7ra,) 2 2—/ n ( ai _ a .\ J 

\ rn v mj v J ' N 

(23) 


If we use the definition in Eq.(7) we get (Cz)jv+i = 
{Ci)n i thereby proving that Ci is an integral of motion. 


C. Floquet Hamiltonian and quasienergy 
eigenstates 


This satisfies C+u = Miip u . By writing H F ip u = u 
and using Eq. (24), we see that the eigenvalue equation 
is satisfied when, w = 27ra • M(mod 27 t). 

We can also compute the momentum-momentum cor¬ 
relator, (piPj) — (Pi)(Pj) for i ^ j over quasienergy eigen¬ 
states. The correlator over this state follows as 


( PiPj > ^ (Pi) (Pj) 


1 V- l^ml 2 

4 ^ 7r 2 (a,; — a,) 2 

m v J ' 


(27) 


In the case of a resonance ai —> aj , this correlator clearly 
diverges, while in the localized case it remains finite. 


VI. EXPERIMENTAL PROPOSAL 

A natural venue for realizing the interacting kicked 
rotor model is superconducting grains. The Hamilto¬ 
nian, Eq. (3) could be implemented using a chain of 
voltage-biased superconducting grains coupled to each 
other using Josephson junctions. Consider a chain of 
grains gated by a ground plane which is resistively con¬ 
nected to ground (see Fig. 2). We then supply a gate 
voltage Vi to each grain, and connect them to each other 
by Josephson junctions J^. 

Because the voltage on the i-th grain is locked to be 
V, its phase winds with an angular velocity (pi = 2 eVi/h. 
Therefore, the resistance and gate capacitor can be ig¬ 
nored when writing the effective stationary part of the 
Hamiltonian: 

I hi = E <1+ - ^2 + C05 #i - + )■ (28) 

i ij 


A special case of our model is when there are no reso- In addition, the “kick” term is obtained by connect- 

nances, namely, all IOMs associated with the momentum ing each grain to a common macroscopic superconductor 
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FIG. 2. Schematic of chain of voltage-biased (Vi) supercon¬ 
ducting grains coupled to each other (J,) using Josephson 
junctions. Each grain is connected ( K ,) to a grounded com¬ 
mon macroscopic superconductor stroboscopically through a 
strong Josephson coupling. 


(which is itself grounded), for a short time and through 
a strong Josephson coupling. 

n*) = -£ K(t — nT) cos (j>i, (29) 

i, n 

with K(t) = K when |f| < r. To make the kick term 
as close to a delta function as possible, we must have 
2ctV.iT <C H2ir, and Kt ~ H, and r -C T. A diagram of 
the circuit for a nearest-neighbor interaction is shown in 
Fig. 2. Identifying, fa = 27ra* = 2 eVi/h and pt = 

hqi/ 2e we see that we indeed obtain the Hamiltonian of 
Eq. (3). 


VII. CONCLUSIONS 

In this work, we introduced the concept of dynam¬ 
ical many-body localization and presented an exactly- 
solvable model of driven linear rotors, which exhibits this 
phenomenon. Although the model possesses a full set of 
integrals of motion, it is shown that dynamical MBL is 
accompanied by the emergence of additional integrals of 
motion, local in momentum space. We believe that this 
observation has important general implications for under¬ 
standing dynamics of interacting many-body systems. 

We have shown that these integrals of motion break 
down due to two types of resonances indicating delo¬ 
calization in momentum space. One type of resonance 
originates from commensuration of the external driving 
period with the parameters of the system and the other 
from the static interactions. An interesting feature of 
this model is that the total energy in the system, that 
is linear in momenta, fails to be a good indicator of dy¬ 
namical localization, since when momentum conserving 
interactions delocalize momentum, momenta of interact¬ 
ing pair grow in opposite directions. 

Moreover, we have shown, (see Appendix) by utilizing 
the lattice mapping introduced by Fishman and Grempel 


and Prange, that our model maps into a disordered lat¬ 
tice with as many dimensions as there are rotors. Based 
on this observation, we argue that what is observed is an 
Anderson type localization, particularly of the type seen 
in correlated disorder systems. 

We also have proposed an experimental setup com¬ 
posed of Josephson junctions and superconducting grains 
to realize the model Hamiltonian. 

Finally, we emphasize that the results presented here 
can apply to more generic non-integrable systems. For 
example, a recent work ’ 2 considers interacting kicked 
Dirac particles with individual Hamiltonians, Hq = 
2tc aa x p + Ma z , and provides a simple argument that 
this non- integrable system also exhibits MBL. First, this 
model also exhibits localization when a’s are generic dis¬ 
tinct irrationals. Second, although the number of inter¬ 
acting particles that can be considered numerically is lim¬ 
ited, note that at large momenta the Dirac model crosses 
over to the linear model considered here. This suggests 
that MBL should be robust to a class of non-integrable 
generalizations, for any number of interacting rotors. 
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Appendix A: Mapping between Quantum Kicked 
Rotor and the tight-binding model 


In the first part of this section, we derive the known lat¬ 
tice mapping of the one dimensional kicked rotor 1 ’. Once 
we establish this derivation, we show that the many-body 
linear kicked rotor (3) also admits a lattice mapping of 
a particle on a d-dimensional lattice. We emphasize that 
existence of such a mapping in the many-body case is 
limited to the case of linear p model (1 = 1). 


1 . Lattice model of single quantum kicked rotor 


In the introduction we considered the time dependent 
Schrodinger equation for a kicked rotor. The kinetic part 
was considered to be p l . For the quadratic kicked rotor 
case l = 2 and the linear kicked rotor that we consider in 
this work is l = 1. Also h = 1 and kicking period T = 1 






























are used in this equation. 

i d t ip(0,t) = [27ra(— idg) 1 + K(0) ^ S(t — n)}ip(d, t). 

n 

(Al) 


The above equation can be solved for if>(9,t) = 

e~ lwt u(Q,t), where the function u has the unit period¬ 
icity of the driving. Let u ± defines the state just before 
and after the kick and are connected by the evolution 
operator in the following way, 

U+ = e - iK{9) u-, u~ = e i ( < *’- 2 ’ r “P , ) tl + (A2) 

Define the following: u = u+ + u ; exp(— \K(0)) = 
(1 — iIF(d)) _1 (l + i W(0)), exp(—i[27rap ; — w]) = (1 — 
iT(0)) _1 (l + i T(d)). Based on the above definitions, 
u^ 1 = (1 =F i T{0))u and 

[T(0) + W(8)]u = 0 (A3) 


is obtained. 

Fourier transforming the above expression, we get the 
following tight binding model, 


^ ^ bF m _ n U n -f- — Eu m . (A4) 

n^m 


Here, the energies and hoppings are: 


1'Em—n 


1 

2-7T 


0 -i(m-n)G 



T m = tan 



2'Kam 



(A5a) 

d<9, 


(A5b) 


This completes the derivation for the lattice mapping 
for the single rotor case. In the following section, we 
generalize this derivation to demonstrate the existence 
of a lattice mapping for the interacting rotor model of 
Eq. (3). 


2. d-dimensional lattice model 


eqs. (8) to (13). The time dependent Hamiltonian in 
Eq. (3) produces the same Floquet operator as, 

d oo 

H ( t ) = 27r otiPi + V F 5{t — n), (A6) 

i—1 n——oo 

where we have defined, 

= E *$) + \ E 4 - § o) (A7) 

*=1 

= E E k ^ im§j + \ E E 'v/d'"' 1 "- "'*• (as) 
2=1 m m 

The factorization enables to treat on site and interaction 
potentials on equal footing. The is defined in Eq. (9) . 
Moreover, we write: 

V F =Y,V?e™- 9 . (A9) 

rh 

The equivalence between Eq. (A7) and Eq. (A9) is con¬ 
ceptually very simple, but somewhat harder to put into 
words. It is probably best to give an example. Take the 
two rotor case i.e. d = 2. Now to are two dimensional 
vectors with integer components. Let us fix m £ Z. If 
to = (m, 0) then V? = k m . If to = (0, to) then Vj[ = k m . 
If to = (to, —to) then Vj[ = b^/2. For all other vectors, 
Vj[ = 0. The summation in Eq. (A9) is over all possible 
vectors to. Generalizing this notation to d dimensions, 
it is possible for both the differences of angles and an¬ 
gles themselves to appear as rh ■ 6. Treating a’s and p’s 
as d-dimensional vectors as well, the Hamiltonian can be 
succinctly as: 


„ OO 

H(t) = 2nafi+^2Vrf l e i ™- 9 ^ S(t-n). (A10) 

rh n =—oo 

Following Ref. [17], the above driven Hamiltonian can 
be mapped on to a d-dimensional lattice model, which 
is closely related to the lattice mapping outlined in the 
previous section for the single rotor case. 

Hfri 1 n'U J n = d^rh^rh T ^ ' dViftftUft = EUfft. (All) 

h 


In this section we show that there exists a d— dimen¬ 
sional lattice model corresponding to the d particle inter¬ 
acting version of the kicked rotor model in Eq. (3) . Such 
a mapping has been previously identified for the case of a 
d rotors driven by an interaction potential 1 '. Notice that 
in our case, the interactions are encoded in the stating 
Hamiltonian H 0 . However, we show that for the linear 
momentum dependence of the kinetic term, static inter¬ 
actions can expressed as the driven interactions and the 
rest of the lattice mapping simple follows from Ref. [17]. 
In order to establish this mapping, we use the factoriza¬ 
tion of the Floquet operator discussed in Sec. (V), see 


Here, to and n are vectors that contain integers that cor¬ 
respond to the quantized eigenvalues of the angular mo¬ 
mentum operator. The hopping and onsite terms are 
defined as, 


Wrri.— n — ESft 


1 

(27r) d 




= tan ( - [w — 27ra • to] 


| d(T 

(Al2a) 

(Al2b) 






Here, we defined E = —(1/[27t] d ) J tan (v F (6)/2^ so as 


to make Wq = 0. 
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